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Abstract We present here results from our search for X-ray pulsations of the neutron star in 
the low mass X-ray binary EXO 0748-676 at a frequency near a burst-oscillation frequency 
of 44.7 Hz. Using the observations made with the Proportional Counter Array on board the 
Rossi X-ray Timing Explorer we did not find any pulsations in the frequency band of 44.4 Hz 
to 45.0 Hz and obtained a 3(t upper limit of 0.47 % on the pulsed fraction for any possible 
underlying pulsation in this frequency band. We also discuss the importance of EXO 0748- 
676, as a promising source for detection of Gravitational Waves. 
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1 INTRODUCTION 

Gravitational Wave (GW) emission is an alluring phenomena, whose detection has eluded the physi- 
cists/astrophysicists so far Direct detection of GWs is technically a very challenging task. Accreting neutron 
stars in Low Mass X-r ay Binary (LMXB) systems, by virtue of their rapid oscillations, are potential sources 



for detection of GWs ( Watts et al. 



2002; 



Haskell et al 



20081). Due t o various reasons, such as, magnetic deformation 



20081) and crustal mountains (IBildsten , 



1998 



Haskell et al 



2006; 



Cutlet 



Vigelius & Melatos 



20081) . these compact objects can develop a quadrupolar asymmetry which leads to GW emission. The GW 



emis sion is also believed to li mit the spin frequency up to which the neutron stars can be spun-up by accre- 



tion (iChakrabarty et al 



. 120031) . In addition to being an important class of potential sources for GW detection, 
any detection of GWs from these sources may also provide important information about the structure and 
geometry of neutron stars. It can also be an important tool to understand the interaction of the neutron star 
with its accretion disk. 

Detection of GWs associated with spinning neutron stars will be an arduous job as it requires a huge 
computational power The major hurdles include, low accretion rate s and uncertainty i n the measurement 



2008 



Watts & Krishnan . 



of spi n and orbital parameters, of most of the astrophysical sources (iWatts et al 
20091). There are sources like Sco X-1, which are quite bright, but the spin and orbital parameters are poorly 
constrained. In case of sources emanating weak signals, the data ought to be folded in order to increase 
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signal to noise ratio. But this folding requires a priori and well constrained ephemeris for spin and orbital 
parameters. If the pulse and orbit ephemeris is accurately known, it will require only a single trial to search 
for gravitational waves. However, when parameters are poorly constrained, a large number of trials are 
required and hence the searches become computationally untenable. 

Among the galactic sources, the LMX Bs harboring rapidly spinning neutron stars are the best target s 



for GW searches with LIGO and VIRGO (lAbadie et al 



201C; 



Abbott et al 



2010; 



Sengupta et al. 



201C). 



Among all LMXBs, the spin and orbital parameters and their evolution are best constrained in the accret- 
ing millisecond pulsars. Hence, least amount of computational resources will be required for the search of 
gravitational waves. But the accreting millisecond pulsars are transient systems, active only during short 
outbursts and the long term mass accretion rate is not high enough for detectable GW emission. The magni- 
tude of long term average flux of most of the accretion powered millisecond pulsars is of the order of 10^^^ 
erg cm~^ s~^, which is much smaller than that of Sco X-1 (of the order of 10"^ erg cm~^ s~^). 



Next best sources ar e the burst oscillation sources and the kHz QPO sources (iWatts et al 



2008 



Watts & Krishnan , 



2009b . But the number of trial searches required in these sources can go up to 35 or- 



ders of magnitude, depending on the number of constrained parameters. 

Among the persistent LMXBs, EXO 0748-6 76 is a promising so urce for detection of GWs. It has been 



studied extensively since its discovery in 1985 (IParmar et al. 



ties of this source, during the operation of Advanced LIGO (ISn 



1985 ). Future monitoring of 



possible with the large area X-ray dete ctors of ASTROS A T jPau] , 



ith et al 



2009; 



timing proper- 



Harrv et al 



2010) wiUbe 



a binary orbit with a period of 3.82 hr (IParmar et al 



20091) . Being an eclipsing source in 



19861) . the orbital period, period evoluti on timescale 



and m id-eclipse time are measured with good accura cy over a time base of more than 20 year s (iWolff et al 



20091) . The sour ce also shows irregular X-ray dips 



the neutron star ( Cottam et al 



companion (IParmar et al 



2002; 



Pearson et al 



(Parmaretal 



1986; 



20061 Oze] . 20061 



Wolff et al 



Munoz-Darias et al 



2002). The mass of 



200% and its 



19861) are also known with reasonable accuracy. 



Since its discovery, EXO 0748-676 has been regularly monitored by various satellites. Most of the 
observations reveal a persistent luminosity of ^ ^ 036-37 (p/7 4 kpc)2 erg s~^. The long term 2-20 keV 



average flux is 2 x 10 erg cm ^ s ^ (IWolff et al 



2008a ). It entered a q uiescent phase in 2008 , with 



a decline of up to 2 orders of magnitude in the flux level 



Wolff et al 



source has also re peatedly shown several Type-I X-ray bursts (IGottwald et al 



2008allb 



Wolff etal 



Hvnes & Jonesl 



198 
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ViUarreal & Strohmaveil (120041) 



Gallowav et al. 



2008). The 



1986 



They 



(l2010h 



20021) . Burst oscillations at 44.7 Hz were discovered by 
associated these oscillations with the spin frequency of the neutron star. However, 
recently found another millisecond oscillations feature at 552 Hz during some of the thermonuclear X-ray 
bursts. The actual spin period of the neutron star is not known. 

In the present work, we have tried to verify whether the 44.7 Hz burst oscillations, seen during the 
thermonuclear bursts, can be associated with the spin frequency of the neutron star. A positive result would 
also allow measurement of the semi-amplitude of the binary orbit, and a monitoring of the spin and orbital 
parameters of this source would make deep search for GW emission possible by allowing a long integration 
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Orbital Phase 

Fig. 1 Folded light curve of EXO 0748—676, obtained from observations made with RXTE- 
PCA. The light curve was folded into 64 phasebins with a period of 13766.78824 s. Two cycles 
are displayed for clarity. 

2 OBSERVATIONS AND ANALYSIS 

The data used in thi s work was obtained from observations made with the Pr oportional Counter Array (PCA: 



Jahoda et al. 



Bradt et al 



d 1993b ), on August 



(Il996h ) on board the Rossi X-Ray Timing Explorer (RXTE: 
17, 1997 (Observation ID: 20082-01-02-00). The total exposure time was 37ks, covering about 3 binary 
orbits. All the five PCUs were ON during the observation. The data was collected from the Science event 
mode (E_125/^s_64M_0_ls), having a time resolution of 125 /^s and 64 binned energy channels. The 
raw events were filtered using the SEFILTER and SEEXTRCT tool of the f tools (version V6.5.1) package of 
the HEASOFT. The average source and background count rate during the observation was '^155 and ~^45 
counts/s, respectively, i.e. the source was about 4 times more intense as compared to the background. The 
fight curve was barycentered using the /too/-FXBARY. 

Figure 1 shows the barycentered light curve of EXO 0748-676, folded into 64 orbital phasebins with a 
period of 13766.78824 s (Wolff et al. 2009). The folded light curve was normalized by dividing the count 
rate by the average source intensity. The folded light curve shows irregular dips along with a total eclipse 
lasting for about 0.04 orbital phase, i.e. 550 s. 

For a neutron star in a binary system, the arrival times of the pulses vary with the orbital motion of the 
pulsar When the pulsar is at Tt,/2 (the farthest end on the orbit from the observer, also referred to as the 
superior conjunction), the photon takes a longer time (by ' ) to arrive (here, i is the source inclination 
angle and siiii is the projected semi major axis). Therefore, assuming a circular orbit for this LMXB, 
the time delay of photon arrival time caused by the orbital motion can be expressed as: 



''emission 



Qx Sim I 27r(t 

H COS 



emission 



orbital 



(1) 
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In order to search for pulsations in EXO 0748-676, we first corrected the arrival time of the photons for 
different assumed values of ax s'mi. Based on the known values of the mass function and the mass of the 



neutron star and its companion jMunoz-Darias et all 120091) . the value of sin i is expected to be around 



730 It-ms. The light curve was corrected for a^. sin i over a wide range of 100-1200 It-ms, with a resolution 
of 2 It-ms which is about 10% of the pulse period around which the search was made. 

After each correction of the light curve for the orbital motion, we searched for pulsations using the 
pulse folding and maximization /too^-EFSEARCH. This is a well established technique and has been 
successfully used to determine the spin period and orbital parameters in various sources (Jain, Dutta & 
Paul, 2007; Jain, Paul & Dutta, 2010). The burst oscillations from EXO 0748-676 have been reported to 
be centered at 44.7±0.06 Hz (Villarreal & Strohmayer 2004). This indicates that the neutron star spin 
frequency probably lies between 0.02234 and 0.02240 s. Using the /too/-EFSEARCH, the light curves were 
folded into 16 phasebins. We searched for 300000 trial periods, with a resolution of 1 x lO^*' s and centered 
at 0.02237 s. This enabled us to search for period over a range which is 5 times larger than the expected 
range for the neutron star spin period with a maximum pulse phase smearing of 10%. Each folded Ught 
curve was fitted with a constant and determined. If the trial period and the a^- sin i value are not correct, 
then the folded profile is smeared and the of the fit will be small, close to the number of phasebins. But 
if the trial period is correct, the pulse profile is reproduced correctly and a constant fit is expected to give a 
high x^- Therefore, the trial period corresponding to a large will represent the true pulse period in the 
light curve, if any. In this process, the a^ sin i will also be determined. 

Figure 2 shows a few samples of the distribution with the trial spin periods for different values 
of ax sini. The value of ax s'mi used for orbital correction is mentioned in each case. We did not find a 
significantly high x^ for any spin period in the entire range of a^. sini. This is shown in Figure 3, where 
a histogram has been plotted which shows the number of times a specific x^ occurs for all the 300,000 
trial spin periods. Histograms were created for the entire range of ax sini and were overlaid. A distribution 
centered at a x^ of about 14-15 is seen in each case. A large deviation from this curve will indicate detection 
of a spin period for the neutron star As seen in Figure 3, a of about 70 is seen in a few cases, which is an 
offset from the x^ distribution curve. But detection significance of these deviations is poor as can be seen 
in Figure 2. 

In order to put an upper limit to the pulse amplitude present, if any, we corrected the light curve for 
the orbital motion with an arbitrary ax sin i. The light curve was then folded with a period of 0.00237 s. A 
model consisting of a sine curve was fit to the folded profile. This gave a 3 ct upper limit of 0.47 % for the 
pulse amplitude. 



3 DISCUSSION 

We have presented results from an X-ray timing analysis of the eclipsing LMXB, EXO 0748-676 in which 
we searched for pulsations at frequencies around the 44.7 Hz burst oscillations in this source observed 
during some Type-I X-ray bursts (Villarreal & Strohmayer 2004). We conclude absence of any pulsations 
around the reported burst oscillation frequency. We have determined a 3(7 upper limit of 0.47 % on the pulse 
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Fig. 2 A sample of variation with trial spin periods for different values of a^; sin i for EXO 
0748-676. 



Recently, Galloway et al. (2010) reported another burst oscillation feature at 552 Hz, which has higher 
detection significance and could be the true oscillation frequency. In order to detect pulsations around 552 
Hz during non-burst episodes, a similar analysis, as explained above is ongoing. 

As mentioned before, detection of GWs is technically, an enormous task. The foremost sources for the 
detection of GWs are the neutron stars of/in LMXBs, where asymmetries are induced in the neutron star 
crust by accretion. The prerequisites for a potential target for GW searches, are, the strength of the signal 
and well constrained orbital and spin parameters. Sources like Sco X-1 and GX 17h-2 may have a large 
GW amplitude, owing to their high X-ray flu x and mass accretion rates 3 x 10^^ Mq yr^^ - 2 x 
10~^ Mq yr~^) (iStrohmayer & BildstenL 120051) . But the orbital and spin parameters are poorly constrained. 
The accretion powered millisecond pulsars, have well constrained orbital and spin parameters, but these 
are tra nsient sources, and have low lo ng-term average accretion rates (with magnitude less than 10^^^ Mq 
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Fig. 3 Histogram showing the variation of for 300,000 trial spin periods of EXO 0748-676 
and for the entire range of ax sin i. 

EXO 0748-676 is definitely a promising source for detecti on of GWs. The l ong-term mass accretion 



rate of EXO 0748-676 is quite high (-2.9 x lO^^" Mq yr'^ dWolffetal 



2005b ) and most of the orbital 



parameters are well constrained. If pulsations are detected, it will make this source an ideal candidate for 
search of GW emission. 
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